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We introduce an analytical iterative method, the density matrix recursion method, to generate 
arbitrary reduced density matrices, and bipartite as well as multipartite entanglement of multi- 
legged periodic spin-1/2 quantum Heisenberg ladders with an arbitrary number of legs. We apply 
this technique to distinguish between even- and odd-legged ladders. Specifically, we show that while 
multipartite entanglement decreases with increasing system-size for even-legged ladders, it does the 
opposite for odd-legged ones. 



I. INTRODUCTION 

The Heisenberg ladder is an interesting platform to in- 
vestigate quantum many body systems in the intermedi- 
ate sector between the one- and two-dimensional lattice 
structures [1]. The possibility of relating doped even- 
legged Heisenberg ladders to high-temperature supercon- 
ductivity [2-4] make such quantum systems extremely 
important. The characteristic pseudo 2-D structure of 
ladders evokes considerable interest in several other areas 
of condensed matter physics [5, G] and quantum informa- 
tion [7]. On the other hand, concepts like entanglement 
and quantum correlations have been applied to under- 
stand critical quantum phenomena in spin systems [n], 
including in Heisenberg ladders [7]. Properties of such 
systems have also been tested experimentally in several 
systems or proposals thereof have been presented, includ- 
ing in compounds and cold gas [!)]. 

A striking feature of the Heisenberg ladder is that 
the interpolation from the 1-D spin chain to the 2-D 
square lattice by gradually increasing the number of legs 
is not straighforward. The quantum characteristics of 
the Heisenberg ladder ensures that the odd- and the 
even-legged ladder ground states have different corre- 
lation properties. "Even ladders" have a finite-gapped 
ground state excitation and exponential decay of two-site 
correlations, while "odd ladders" are gapless and have 
power-law decay [G, 10]. The results of Heisenberg lad- 
ders cannot therefore be directly extrapolated to the 2-D 
regime. This difference in quantum characteristic of odd- 
and even-legged ladders may lead to interesting features 
in the entanglement properties of the system. However, 
calculating the bipartite or multipartite entanglement in 
large-sized multi-legged Heisenberg ladder ground states 
is numerically challenging. Within a dimer-covering ap- 
proach [G], it is possible to derive energy density and 
spin-correlation functions for two-legged ladders by using 
generating functions [11] or state iterations [12]. Approx- 
imate solutions may also be obtained for the four-legged 
ladder [13]. 

In this paper, we introduce an analytical iterative tech- 
nique, the "density matrix recursion method" (DMRM), 
to obtain the reduced density matrix of an arbitrary num- 
ber of sites of a quantum spin-1/2 Heisenberg ladder with 



an arbitrary number of legs and with both open and pe- 
riodic boundary conditions. Specifically, within a dimer- 
covering approach, we find separate iterative formalisms 
for even and odd ladders. These partial density matri- 
ces can be used to calculate the scaling and behavior 
of single-, two-, and multi-site physical properties of the 
whole system, including two-site correlations and bipar- 
tite as well as multipartite entanglemeirt. We then apply 
our method to obtain the nature of a genuine multipar- 
tite entanglement, quantified by the generalized geomet- 
ric measure (GGM) [14]. Multisite entanglement is a use- 
ful tool in the topological study of ground states of spin 
systems [15]. Using a multisite entanglement measure is 
advantageous over the bipartite correlation measures, as 
the former involves a greater distribution of information 
between the multisite sub-regions and the hidden corre- 
lations due to tracing out is not negated. We compare 
odd- and even-legged ladders by using genuine multisite 
entanglement, and find that the GGM of odd-legged lad- 
ders increases with system-size while it decreases in the 
even ones. 

The paper is organized as follows. We introduce the 
model state in Sec. H. The genuine multiparty entangle- 
ment measure is defined in Sec. HI. The main results are 
presented in sections IV and V. In particular, the density 
matrix recursion methods for even and odd ladders are 
respectively introduced in sections IV A and IV B. We 
present a conclusion in Sec. VI. 



II. MODEL STATE 

Within a dimer-covering approach, the ground state of 
a antiferromagnetic Heisenberg ladder consists of near- 
est neighbor dimer pairs on the spin-1/2 lattice. Each 
bipartite lattice site is occupied by a spin-1/2 qubit. 
The pairs are limited to nearest neighbor (NN) directed 
dimer pairs, as governed by the interactions on a bipar- 
tite lattice. The multi-legged Heisenberg ladder consists 
of coupled parallel Heisenberg chains with M sites on 
each chain labelled from 1 to Af. The number of chains, 
referred as legs, of the ladder are labelled as 1 to M. 
The total number of spin-1/2 sites in the entire bipartite 
lattice is s (s = MM). The vertical chains, referred as 
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rungs, are labelled from 1 to n, similar to the sites on 
the chain, and each rung contains M spins. See Fig. 1. 
The periodic ladder is conditionally defined by allowing 
dimer states to be formed between rungs 1 and n [11]. 
Of course as numbers, n = Af. The quantum state un- 
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der study is therefore [17] |7V) = J2k 

where \{ai,bj)) refers to a dimer between sites at and bj 
on the bipartite lattice. are NN sites with i ^ j- A 

product of all such dimers on the spin-1/2 lattice consti- 
tutes a complete dimer covering. The summation refers 
to the superposition of all such dimer product states or 
coverings. 
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FIG. 1. (Color online) The M-legged, A/'-runged Heisenberg 
ladder in the form of a bipartite lattice. The solid circles are of 
sublattice A, while the hollow ones are of B. The red arrows 
(solid) show the nearest-neighbor dimer states {\{ai,hj))) from 
a site in sublattice A to another in B. The arrows (dashed) at 
the boundary indicate that the lattice is with periodic bound- 
ary condition. The figure also shows the two-site ladder state, 
|2), and the one-site rung state, jl), which will be profusely 
used in the iteration method presented. 



III. GENUINE MULTIPARTITE 
ENTANGLEMENT 

To quantify the amount of genuine multisite entangle- 
ment present on multi-legged ladders of different sizes, 
we consider a genuine multipartite entanglement measure 
called generalized geometric measure (GGM) [14, 18]. 
The GGM of an iV-party pure quantum state |c/)7v) is 
defined as 



1-A 



where Amaxd'/'Af)) = max |(x|<?!>Ar)|, with the maximiza- 
tion being over all pure states \x) that are not gen- 
uinely A^-party entangled. An A^-party pure quantum 
state is said to be genuinely multiparty entangled if it 
is entangled across every bipartite split of the N par- 
tics. It was shown in Ref. [14, 18] that £{\4)m)) = 
1 - iaax{X\.g\A[J B = {1, 2, . . . , A^}, A n B = 0}, where 



\a:B is the maximal Schmidt coefficients in all possible 
bipartite splits A : B oi [^at). The GGM can therefore 
be calculated once the reduced density matrices of the 
ladder states can be calculated. 



IV. DENSITY MATRIX RECURSION METHOD 

The iterative method for generating arbitrary local 
density matrices are introduced in the following two sub- 
sections for the even and the odd ladders. We will show in 
a subsequent publication that both the iterations devel- 
oped here can finally lead to obtain an iteration method 
for the 2D square lattice [19]. 

A. The Even Ladder 

The recursion relation for the non-periodic (i.e. with 
open boundary conditions) even-legged ladder can be 
written as 

\N + 2) = |A/-+l)|l)„+2 + |AA)|2) 

= |A/')|2)„+i,„+2 + lA/"- l)|2)„,„+i|l)„+2,(l) 

where |7V) is an even-legged ladder containing N sites 
in each chain and corresponds to the rungs numbered 
from 1 to n. On the other hand, e.g. \M)2,n+i repre- 
sents an even-legged ladder containing N sites in each 
chain, and corresponds to the rungs numbered from 2 
to n 4- 1. |1) is a vertical chain containing M spins 
representing a single rung of the ladder with the dimer 
coverings between adjacent sites. |2)„-|.i^„-|_2 is a two- 
rung even-legged ladder of rungs n 4- 1 and n 4- 2. Also, 
|2)„+i,„+2 = |2)„+i,„+2 - |l)„+i|l)„+2. If no site num- 
bers are mentioned in subscript of [A/"), it is to be con- 
sidered to be from 1 to n = M . 

The periodic boundary condition entails that the lad- 
der also form dimer states between the rungs 1 and n. 
Accounting for the additional states in the ladder sys- 
tem, due to the periodicity, the recursion relation can be 
modified into 



\M+2)p = |A/' + 2)i,„+2 + |A') 



+2,1, 



(2) 



where the subscript P stands for a periodic ladder state. 
Throughout the paper, a state without a subscript P will 
imply a non-periodic state. The density matrix for the 
periodic ladder state can be calculated by using (2): 



AN+2) 
Pp 



\M+2){Af+2\p 
- |A^ + 2)(A^ + 2| + |AA + 2) (AA|2,„+i(2|„+2,i 
+ |AA)2,„+i|2)i,„+2(A/'+2| 
-H|A/-)(A/-|(2,„+i)|2)(2|(i^„+2) (3) 

If, for example, we trace out all but two rungs (say n 4- 1 
and n 4- 2) from the state, we will obtain a two-rung 
(mixed state) ladder containing 2M spins, where M is the 
number of legs of the ladder. Note that M is even here. 
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These reduced states can then be used to calculate the 
multisite entanglement properties of the ladder system. 
In particular, 

/'pd+i.n+a) =tri..„(|AA + 2)(AA + 2|p) (4) 

is the two-rung state of the periodic ladder. Now, 
using relation (1), we can obtain the trace for the 
non-periodic part to have tri n{\J^ + 2) (A/" + 2|) = 
tri^.„[[A/-)(A/'||2)(2|(„+i,„+2) + W - 1)(AA - 
l||2)(2|(„,„+i)|l)(l|(„+2) n+l.n+2 
l|(2|„,„+i(l|„+2 + lA/"- 1)|2) 

n,n+l I ^ / n-\-2 

Tracing over the rungs 1 to n, we get 

(2) 

Pin+l,n+2) = ^A/'|2)(2|(„+l,„+2) + ^AT-lPn+l ® |l)(l|(„+2) 
{l\n+2{^Mn+l + h.C.) (5) 

where Zjv = (AAIA/"), Pn+i = tr„[|2)(2|(„_„+i)], and 
{^Af\n+i = (2|n,n-i-i(A/' — l\Af). Extending the trace to 
the periodic ladder state, we get 

Ppl+l,n+2) = Pnil,«+2 +tri..„[|AA)(AA|(2,„+i)|2)(2|(i,„+2) 

+ (|AA)2,„+i|2)i,„+2(AA+2| + /i.c.)] 

= Pnll,n+2 + /3l(ji+l,n+2) + iP2{n+lji+2) + ^'C-) (6) 

Here, /3^(„_^i_„_^2) = 2aA-i|1)(1|(„+i) ® Pn+2 + 

Z^_2P7i+l O Pn+2 + _(|l)(?A^-l|(n+l) ® Pn+2 + h.c), 
with ($AA-l|n+l = (2|,i,„+l(A/' - 2|2..n-lA/' - 1)2. .n- 
And, f^2{n+l,n+2) ^ 1 2) „+l,„+2 (1 |n-|-l (^A^' |ri+2 + 
\'2')n+l,n+2{<t'f\f\n+l,n+2 + Pn+1 ^ I l)n-|-2 ('Ca/'-I |n+2 + 
3lCl)n+lJl)n+2(l|«+l(W-l|n+2, whcrC (0Ar|„+l,„+2 = 
(2|i,„+2(2|„.„+l(7V - 2|2.„-l7V), {m-i\n+2 _ 

(2|i,„+2(A/'-l|2,„AA-l)i,„-i|l)„,and^= - 

The local density matrices can be calculated by using 
the iterations, derived below, of the complete and partial 
inner products defined above. Using (1), the normaliza- 
tion of the non-periodic ladder is given by 

= (AA|A/-) 

^AZj^^i+BZ^^2 + 2Cy]^_2 + 2Vy^_^. (7) 

where_^ = and iS = (2[2)_. C, V, C, and_2? arc given 
by {1|2) :=C|1)+P|1) and (1|2) =C|1)+P|1), where |1) 
is a vertical rung with a periodic dimcr covering between 
the topmost and lowermost sites of that rung. The other 
terms in the normalization can be calculated as follows: 

y)^ = (AA|AA- = Az^^i + cy'^_, + vy^^_,, 
ylr = - = Az^^i + cy]^_, + vyl_,. 

The other term in the expression for /0„_j.]^ ^+2 that is to 
be calculated using iterations is 

{iM\ = {%n.n+l{M-l\M) 

= {l\{CA)^ + CA^j^) + {l\{VA\ + 'DA^j^), (8) 

where the iteration variables are given by 
TV M 
^Jv = X! ^A/- = X! Z^^-^h^_l, (9) 

i=l i=2 



with = Cg^+Cht, hi+i = Vgi+Vhi, go = 1, /ig = 0. 
The extra iterative terms in the periodic two-rung dcn- 

(2) 

sity matrix, Pp^^^^i n+2}^ '^^^^ ^® expressed as 

{4'^/\n+l,n+2 — (2|l,n+2(2|n,n+l(A/' — 2|2,n-lA/') 
= (2 1 1^)1+2 (2 In, n+1 (-^1^11)1 1 l)n 
+ A'2^|i)l|l)n + A'3^|l>l|l)n + |l)l|l>n 
+ (2|n-2,n-l...(223|2>12...|2)n-l,n). (10) 

Here, the iterative variables, Xf^ , are defined by using 



(9), in 


the following way: 










= ^92i{A]^^l^2i + 


' C^7V-2-2j) " 


h g2iVA\f_2_ 


-2'i> 


x^ = 


= ^92i{Alf^^_2i 4 


' CA^_2_2i) - 




-2'iJ 


xf^ 


:^/l2,(A^-l-2z^ 


- CA\_2_2i) ' 


f h2iVAlf_2. 


-2ii 


yN 
'T-4 — 


= ^/^2^(A2^_1_2z^ 


- CA\f_2_2i) - 


f h2iVAl^_2_ 


-2i! 



(11) 



where the summation is from i = to A/". The variables 
5i, hi mimics the same relations as gi,hi, with initial con- 
ditions 5o = 0, /iQ = 1. Similarly, we can derive 

(W-l|«+2 = (2|l,„+2(A/'- Ib.nA/"- l)l,n-l|l)n 

= {AA'^_^ + ^A^„i)(C(l| +P(1|) + (2|„+24 

(^-1^-1(112(1^, + ;f2'^-'(l|2(l|n + ;f3^-' 
(l|2(l|n+A'f-l(l|2(l|n)|2)l,2|l)n, (12) 

where the iterative variables can be derived using (11). 
These iterations can be used to obtain the partial den- 
sity matrices of a periodic even-legged ladder, and hence 
its bipartite as well as genuine multipartite entanglement 
and other single and multi-site physical quantities, pro- 
vided the values of A, A, B, C, C, V, V, Zi, yl, and yf 
are exactly calculated. 

The size of the reduced density matrices depend on 
the type of the ladder considered. For example, a two- 
rung reduced density matrix for an M-legged ladder is 
a T = 2Af-spin matrix. Other reduced density matrices 
of spins smaller than t, required e.g. for calculating the 
GGM, can be obtained from the r-spin matrix by partial 
tracings. 



B. The Odd Ladder 

The periodic recursion for the odd ladder is rather dif- 
ferent than the even ladder. This is due to the fact that 
there exists no analogous state for |1) in the odd ladder. 
The non-periodic |7V) can be written as a series of A/" = 2 
odd-legged ladders. \J\f) = |2)i,2|2)3,4...|2)„-i,„. The pe- 
riodic recursion for the odd Af-legged A/'-spin-pcr-chain 
ladder is then given by 

l,n|2)n+l,n+2 2,n+l|2)n+2,l (13) 
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where \2)n+i,n+2 is a two- rung odd-lcgged ladder at rungs 
n + l,n + 2. The density matrix can be calculated, as 
before, using (13): p^^^ = \J\f + 2){Af + 2\p. The re- 
duced density matrices can be obtained by tracing out 
the requisite number of spins. In particular, for obtain- 
ing a two-rung reduced density matrix, we trace out the 
rungs ranging from 1 to n: 

Ppln+i.n+2) = tri..„[|A/' + 2)(A/- + 2|p] 

= tri..„[|AA)(AA|i,„|2)(2| n+l,n+2 
+ |A/')(A/'|2,„+i|2)(2|„+2,i 

+ {\-^) l,n\'2) n+l,n+2 

After simplification, the above equation reads 



(2) 

^P(n+l,n+2) 



^A/'|2)(2|(„+i,„+2) + Z^f^2Pn+l ® Pn+2 
(|2)„+l,„+2(f^AAU-|-l,n+2 + h.C.) (14) 



where {ilu\n+i,n+2 = (2|i,„+2(A/'|2,«+i|A/')i,„, and the 
normalization is given by 



N/2 



(15) 



Now, for a Af-legged ladder, Z2 corresponds to a two- 
legged A/-site-per-chain ladder, and can be calculated 
from the previous section. The recursion for {^j\f \ is 



{^N'\n+l,n+2 — (2|l,ri+2(A/'|2,n+l|A/')l,„ 

= (2|l,n+2(2|2,3---(2|n,n+l|2)l,2|2)3,4...|2)„_i^„. 



(16) 



Hence the computation involves writing an algorithm 
to iterate the step (2|i^i_|_3 (21^+1^^4.2 |2)i^i+i. Once the 
reduced density matrices are obtained by the iterative 
method, we can again use them to obtain the different 
single- and multi-site physical quantities of the system. 



V. EVEN VERSUS ODD 

To illustrate the effectiveness of the density matrix re- 
cursion method, we apply it to obtain a multisite en- 
tanglement of multi-legged ladders with both even and 
odd number of legs. In particular, we consider two- and 
four-legged ladders among even ladders, and three- and 
five-legged ladders among odd ones. The iterative vari- 
ables can be evaluated by using their explicitly calculated 
initial set of values. The iterations provide the reduced 
density matrices of the system, which are thereafter uti- 
lized to obtain the GGM. 

Specifically, for M = 2, the relevant initial parameters 
a_re Zq = 1, Zi = 2, ^ = 2, ^ = 2, C = 1, 2? 0, C = 0, 
2? = 0, yl ^ 2, and yl = 0. For M ==_4, the initial 
parameters are Zq — 1, Zi ~ A, A = A, A ~ 2, C = 5, 
V = 1, C = 2, V = 3, yl ^ i, and y^ = 2. A similar 
analysis can also be done for higher even-legged ladders. 

In a similar fashion, the density matrix recursion 
method can be applied to obtain the multisite entangle- 
ment for odd ladders. The initial parameter required in 
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FIG. 2. (Color online) Genuine multisite entanglement de- 
creases with system-size for even ladders. We perform the 
iterations for M — 2 and 4. The iterations are carried out 
until 18 rungs, i.e. 36 and 72 sites respectively for M — 2 and 
4. However in both the cases, the GGM converges much be- 
fore those sizes. The vertical axis represents the GGM, while 
the horizontal one represents the number of rungs. Both axes 
are dimensionless. 




FIG. 3. (Color online) Genuine multisite entanglement in- 
creases with system-size for odd ladders. We perform the 
iterations for M = 3 and 5. The iterations are until 18 rungs, 
i.e. 54 and 90 sites respectively for M = 3 and 5. Again in 
both the cases, the GGM converges much before 18 rungs. 
The vertical axis represents the GGM, while the horizontal 
one represents the number of rungs. Both axes are dimen- 
sionless. 



the recursion for M = 3 is Z2 44, and for M = 5 is 
Z2 = 804. 

We are now ready to compare the multisite entangle- 
ments for the even- and odd-legged ladders. The GGMs 
obtained by the iterative methods clearly capture the 
characteristic complementary nature of even and odd lad- 
ders. We show that the GGM decreases with the increase 
of system-size in the case of even ladders (Fig. 2). The 
opposite is true for odd ladders - the GGM increases with 
system-size (Fig. 3). Although the comparison is made 
by taking M = 3 and 5 among odd ladders, we have actu- 
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ally performed the computations also for M = 1, which 
again shows the characteristic increasing GGM of odd 
ladders. In all the cases, the GGM is calculated by con- 
sidering reduced density matrices upto 2M spins: numer- 
ical exact diagonalizations corroborate that considering 
upto 4 spins is already enough. We have performed the 
iterative algorithms upto 18 rungs in all the cases (which 
e.g. is equivalent to 72 spins for the four-legged and 90 
spins for the five- legged ladder). As seen in the Figs. 2 
and 3, the GGM has already converged much before the 
maximum number of rungs that we have considered. 

VI. CONCLUSION 

We have introduced an analytical iterative technique, 
the density matrix recursion method, which can be effi- 
ciently used to obtain arbitrary reduced density matrices 



of the ground states of spin-1/2 antiferromagnetic quan- 
tum Heisenberg ladders with an arbitrary number of legs, 
within a dimer-covering approach. This technique imme- 
diately allows us to obtain single- and multi-site physi- 
cal properties of the system. In particular, we use the 
method to obtain the scaling of genuine multisite entan- 
glement in both odd- and even-legged ladders. We find 
that the genuine multisite entanglement can capture the 
disparity between even and odd ladders. 
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